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ABSTRACT 

We report the complete photometric results from our Herschel study which 
is the first comprehensive program to search for far-infrared emission from cold 
dust around young brown dwarfs. We surveyed 50 fields containing 51 known 
or suspected brown dwarfs and very low mass stars that have evidence of cir- 
cumstellar disks based on Spitzer photometry and/or spectroscopy. The objects 
with known spectral types range from M3 to M9.5. Four of the candidates were 
subsequently identified as extragalactic objects. Of the remaining 47 we have 
successfully detected 36 at 70/xm and 14 at 160/im with S/N greater than 3, as 
well as several additional possible detections with low S/N. The objects exhibit 
a range of [24]- [70] micron colors suggesting a range in mass and/or structure of 
the outer disk. We present modeling of the spectral energy distributions of the 
sample and discuss trends visible in the data. Using two Monte Carlo radiative 
transfer codes we investigate disk masses and geometry. We find a very wide 
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range in modeled total disk masses from less than 10 -6 M up to 10~ 3 M with 
a median disk mass of order 3 x 10~ 5 M , suggesting that the median ratio of 
disk mass to central object mass may be lower than for T Tauri stars. The disk 
scale heights and flaring angles, however, cover a range consistent with those 
seen around T Tauri stars. The host clouds in which the young brown dwarfs 
and low-mass stars are located span a range in estimated age from ~l-3 Myr 
to ~10 Myr and represent a variety of star-forming environments. No obvious 
dependence on cloud location or age is seen in the disk properties, though the 
statistical significance of this conclusion is not strong. 

Subject headings: protoplanetary disks — stars: formation — brown dwarfs 



1. Introduction 



The existence of circumstellar disks around low-mass stars and substellar objects is now 
a well-established phenomenon. Many of the ob servations of the lo west mass objects have 



taken advantage of the Spitzer Spa ce Telescope ' s (IWerner et al.ll2004j) exquisite sensitivity in 
the 3 - 24/xm spectral region, e.g. lAllers et al.l (120061 ); iHernandez et al.l (120071 ); iGuieu et al. 



( 120071 ): iLuhman et al.l (120051 ) (and many of the other references in Table [3]) to find thermal 
emission above that expected from the photospheres of these objects. With the likely surface 
density gradient of such disks, however, most of the dust is typically located at radii too 
large to emit strongly at those Spitzer wavelengths, particularly for low-luminosity sub- 
stellar objects. At longer wavelengths, 70 - lOOOAtm, few disks have been detected around 
brown dwarfs, e.g. (IKlein et al.ll2003l ; IScholz. Jayawardhana &: Woodll2006t iBouy et al.ll2008 



Riaz fc Gizisl 120081 1 . In addition to being more sensitive to the total disk mass, observations 
at the longer wavelengths are also sensitive to the disk scale height, dust settling, and degree 
of flaring. The relationship between disk mass and structure and planet formation is highly 
uncertain, but nominally one might exp ect objects with lower mass disks to form planetary 
systems with lower average mass planets ( iMordasini et al.ll2012l ). Below some disk mass limit, 
though, there may be too little material for any planet formation, so understanding disks 
around the lowest mass objects will place important constraints on likely planet formation 
around them. 



* Herschel is an ESA space observatory with science instruments provided by European-led Principal 
Investigator consortia and with important participation from NASA. 
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The relatively small number of detections of disks around very low-mass stars and brown 
dwarfs beyond the range of Spitzer's 24/mi MIPS band and its IRS instrument, has made 
it difficult to investigate the properties and evolution of the cool dust component of these 
disks where most of the mass is located. We describe here the results from a Herschel GT1 
(Guaranteed Time, Phase 1) program that obtained 70 and 160/im photometry of a statisti- 
cally significant number of confirmed or candidate brown dwarfs and low-mass stars close to 
the sub- stella r limit . Our early results on the first three objects observed are described by 



Harvey et al.l (120121 ) (hereafter, Paper I). In the following sections we describe the observed 



sample, the details of the observing parameters and data reduction, the overall properties 
of the spectral energy distributions (SED's), and extensive radiative transfer modeling of 
the SED's to extract physical properties of the disks. We conclude with a discussion of the 
overall results from this study based on the complete sample together with existing data in 
the literature on cold dust in young brown dwarf and T Tauri disks. The most important 
goals of this program were to use the Herschel far-ir photometry to place useful limits on the 
outer disk properties around brown dwarfs, in particular, the mass and disk geometry. We 
also aimed at providing as complete as possible SED coverage for future characterization at 
longer wavelengths with the Atacama Large Millimeter Array (ALMA). 



2. The Sample 

Our sample of 50 targets was chosen from the literature to include young objects at or 
below the sub-stellar limit in the nearby regions of star formation listed in Table [TJ Both 
confirmed and candidate BD's were included. Table [1] lists rough ages for each of the clouds 
based on the range of measurements in the literature. Although there are often more precise 
numbers quoted for median cloud ages in the literature, the age spreads observed around the 
median have led us to consider only two age bins in our sample, the clouds with ages < 5 Myr 
and those with ages of order 10 Myr. Because the sub-stellar limit is only roughly defined 
as a function of spectral type and age, we included objects generally with spectral types 
later than M5 with the exception of several in the range M3 - M4.8. The most important 
additional selection criterion was evidence for circumstellar dust emission, likely a disk, as 
shown by 8-24/im excess emission over the photosphere. Our goal was to reach flux limits 
at 70/^m low enough that we would be likely to detect disk emission based on the shorter 
wavelength Spitzer data or to set useful upper limits. Our sample was specifically trimmed 
to not include objects that might be bri ght enough to be d etected in Herschel 's large scale 



shallow surveys of star-forming regions (lAndre et al.ll2010l ). The selection process involved 



"predicting" the 70//m flux density based on the 8-24/zm color and 24/im brightness using 
typical SED's for brighter young T Tauri stars as a guide. As might be expected, the final 
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results showed our rough "predictions" to be accurate to only ± a factor of 3. An additional 
selection criterion was an attempt to select objects that did not lie in regions of strong, 
structured diffuse emission as based on Spitzer 24 and 70/xm images when available. Even 
with this criterion, one object, CFHTWIR96, was completely hidden in much brighter diffuse 
emission at both observed wavelengths, and a number of sources were invisible at 160/im 
against the Galactic background/foreground. 

Because the sample included candidate BD's, a few objects included in the observed 
sample have turned out not to be stellar or sub-stellar objects with disks, including four 
found to be extragalactic and one that may be at an earlier evolutionary state with a 
substantial envelope or extragalactic (ChaJ11080397-7730382). The four extragalactic ob- 
iects include 2MA SSJ16080175-3912316 and IRA CJ16084679-3902074 from the study b) 



lects include 2MA SSJibU«Ui^-dUi^dib and iKA UJibU84b^-dMU2U74 irom tne study by 
Allen et aD (2007 )and found to be extragalactic by Comeron. Spezzi &: Lopez Marti (20091) 



and 



Comeron! (1201 11 1. The source 2MASSJ1129140- 7546256 was found bv lLuhman fc Muench 



(120081 ) to exh ibit emission lines characteristic of a galaxy. The object 2MASSJ16084747- 



3905087 from I Allen et al.l (120071 ) has been found to consist of two nearby objects, one of 
which peaks in the near-ir and the other at longer wavelengths by Luhman (private comm. ) . 
This is most obvious in the Spitzer IRAC band 1 (3.6/im) image where the object is clearly 
double, consistent with it being a blend between a star and extragalactic object. 

Table [3] lists the complete sample of low-mass stars and brown dwarfs, including one 
object, IS0143 that was included serendipitously in the field of IS0138 and which lies at 
the boundary between stellar and sub-stellar objects. With the inclusion of this object and 
the removal of the four extra-galactic sources, the final sample consists of 46 likely or con- 
firmed BD's/low mass stars with circumstellar disk emission, and one object, ChaJl 1080397 



77303 82 (hereafter chal 10804), classified as a "possible Class 1 source" by iLuhman et al. 



(120081 ). We include this latter source because its Spitzer IRAC band 4 (8.0/xm) image shows 
faint nebulosity around the object that is at least consistent with it being a Galactic source. 
We have also measured upper limits to its 1.6 and 2.2/im fluxes from the 2MASS images of 
the area, since it was listed as undetected by the survey Figures [1] and [2] display the physical 
properties of the observed sample graphically. Figure [1] shows a plot of the observed total 
luminosities (1-160/xm) versus spectral type for the 43 objects with published spectral types 
assuming the distances listed in Table [TJ This figure and several later figures show individual 
points with symbols for each host cloud; in addition, the two oldest clouds, Upper Sco and 
TWA, are indicated with red symbols to highlight any possible age-dependent trends using 
the two age bins discussed above. Figure [T] illustrates graphically the wide range in lumi- 
nosities present in our sample. It is also clear that there is a substantial range in luminosity 
at any given spectral type and likewise a range in spectral types at any given luminosity. 
There is no obvious correlation with cloud age for either of these distributions. Figure [2] 
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shows histograms of the observed total luminosities of all the objects and the distribution 
of spectral types for the 43 objects with published types; the separate contributions from 
each host cloud to the total sample are indicated, illustrating the range of luminosity and 
spectral type present in the sample from each host cloud. 

Because of the various features of the above selection process, our study is not useful 
to estimate disk frequencies for young brown dwarfs. The sample, however, is useful to 
compare the frequencies of cold dust emission relative to warmer dust emission due to, for 
example, differences in outer disk mass or structure in the outer disk. More generally, our 
far-ir photometry has the potential to provide the first statistically meaningful sample of 
disk masses around BD's. Since we have objects from various star-forming regions covering 
a range in estimated ages from ~l-3 Myr to ~8-10 Myr, this sample may also be useful for 
investigating disk evolution around BD's, e.g. mass and geometry. We discuss all of these 
issues more fully in £j6l 



3. Observations and Data Reduction 

Our program, GTl_pharve01_2, was designed to ta ke advantage of the excellent sen- 



sitivity and spatial s ampling of the PACS photometer (IPoglitsch et al.l 120101 ) on Herschel 



(IPilbratt et al.ll2010l ). All the observations were made using exactly the same mode with 
essentially the same total integration time. We utilized the PACS "mini-scan-map" mode 
which images simultaneously at 70 and 160/im and uses two separate scan maps at an angle 
of 40° for good 1/f noise reduction. This mode provides high and relatively uniform sen- 
sitivity over an area of order 60"x90" in the center (within 80% of peak coverage). Table 
H] lists the relevant parameters for all the observations reported here. The OBSID's for all 
the observations are listed in Table [3j The la noise levels for the combined OBSID's for 
each object were expected to be ~1 mJy at 70/im and 2.3 mJy at 160/mi according to the 
Herschel observation-planning tool HSPOT which is based on standard aperture photometry 
estimates. In fact, we achieved 10-20% better noise levels at 70/zm probably because of our 
use of psf-fitting photometry instead of aperture photometry. The uncertainties at 160/im 
were almost entirely set by the level of diffuse emission present in the images, which differed 
widely between objects. 

The data were first processed with the Herschel Interactive Processing Environment, 
HIPE version 7.3, with a script that utilizes standard high-pass filtering for point-source 
observations, and which processes the two separate OBSID's for each object together. The 
portions of scan legs where the telescope was stopped or not scanning at a uniform velocity 
were not used for this processing. The output of this script consists of fits files of the image, 
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coverage, and uncertainty. The uncertainty images are not yet reliable at this stage in the 
Herschel mission, so we estimated uncertainties in additional ways described below. 



The next stage of data reduction utilized the ps 



inally developed by the c2d Spitzer Legac y Team (|Harvev et al 



-fitting photometry tool, c2dphot, orig - 



2006 



Evans et all 120071 1 



and based on the earlier DOPHOT tool (ISchechter. Mateo, fc Sahal 119931 ). This tool can 
be used in various modes, including finding peaks above the background and fitting a psf 
to the local maxima, or fitting a psf to any fixed position in the image. This latter mode 
is especially useful for estimating the noise and determining upper limits. The noise levels 
quoted for our sources were determined this way by fitting a grid of nine positions on and 
around the BD position within the high-coverage area. Upper limits were also estimated 
this way as well as by inserting artificial stars into the images, particularly for the 160/im 
data. We also computed aperture fluxes for comparison with the psf-fitted values, and those 
values agree well. In addition, we have computed aperture-flux "curves-of-growth" to be 
sure that there are no systematic problems with the psf-fit photometry. At 160/im most of 
the images were contaminated substantially by structured, diffuse emission that dominated 
the flux uncertainties. 



4. Observational Results 



Figures [3H7] (electronic version only) show images of all 47 of the observed low mass 
young objects. A summary of our measured flux densities is given in Table [31 We have 
not included the four objects fou nd to be extragalactic , but do list the object chal 10804 
described as "possible Class 1" by iLuhman et al.l (120081 ) since its luminosity is probably in 
the same range as the other objects. The entries in Table [3] given as upper limits have been 
determined by finding the maximum brightness source that could have been present at the 
position of the object without us reliably detecting it; therefore, these should be viewed as 
firm limits, equivalent to ~ 3a. For objects where an examination of the image showed a 
possible detection, we report the formal flux level derived from the psf-fit and the uncertainty 
determined as described above using multiple psf-fit positions in blank parts of the image. 

All of the objects in Table [3] have been obser ved over a wide range of wavelengths by 
some o r all of the 2 MASS Jskrutskie et al.l hood ) and DENIS surveys, the recent WISE 
survey (IWright et al.ll2010l ). numerous Spitzer programs using IRAC and MIPS, as well as 
36 with the IRS instrument on Spitzer. This enables us to present very complete SED's for 
all the sources, and as discussed below, enables model fitting to derive limits on physical 
parameters of the circumstellar disks. We show these in two forms. Figure [8] shows the SED's 
of all the objects superposed and normalized to the 1.6/im flux density (or upper limit in 
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the case of chall0804). Figures 191 - [121 show the individual SED's for all the objects except 
chal 10804. 

Figure [8] shows that the observed SED's between 1 and 100/im differ considerably, and 
it suggests also that there appears to be a relatively continuous distribution of slopes in 
the mid- to far-infrared. We can also illustrate this conclusion with the color-color diagram 
in Figure [13] that shows a plot of the ratio of uF u at 24/zm versus 70/mi compared to the 
comparable ratio at 3.6/xm versus 24/xm. This also shows a relatively continuous distribution 
of colors, although 90% of the objects have z/F J/ (24//m) > z/F„(70/mi). The very red object 
chall0804 is clearly beyond the end of the distribution of the other objects that are believed 
to have circumstellar disks but not substantial envelopes. 

One question we can consider with just this simp lied analysis of the SED's is whether 
there is any evidence for differences in disks with age or environment. The symbols for each 
object in Figure [13] have been coded to denote the host cloud in which they are found. Within 
the relatively small number statistics for some of the regions, there is no clear trend in far-ir 
excess as a function of the location/age of the cloud in which each object is found. The 
same conclusions holds if we restrict the sample to a subset of objects with well-determined 
spectral types in the M5-6 range which has a reasonable number of sources. 



5. Radiative Transfer Modeling 



In Paper I we showed results from modeling the disks around the first three objects ob- 
served in this program, 2MASS120733, IS0138, and SSSPM110209. Despite the degeneracy 
inherent in such modeling, we found that for those three cases with excellent SED coverage 
we were able to constrain the circum-object disk masses and geometry to some degree. For 
this current study we performed radiative transfer modeling of the complete set of objects, 
though the results are probably most useful for those with the most complete SED coverage. 
We have not tried extensively to fit the Spitzer IRS data in the area of the sili c ate fe ature 
which is generally viewed as arising in the disk atmospheres, e.g. lOloffson et al.l ( 120091 ). and 
we have assumed a single, uniform grain composition. 



by. 



We used the two radiative transf e r code s discu ssed already in Paper I, MC3 D described 



Wo 



(12006 



f. Henning fc Steckluml (119991 ): WolJ (120031 ) and MCFOST described by IPinte et al. 



20091 ) . Both codes are three-dimensional, r adiative transfer codes using the Monte 



Carlo method. NextGen stellar atmosphere models ( jHauschildt. Allard fc Baron 



1999; Allard fc Hauschildt 



19951 ) were used for the stellar input to the disk models with the photospheric temper- 
ature set by the assumed spectral type where available. With each code we modeled 
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each SED with a parametric disk. We use a dust density distribution with a Gaussian 
vertical profile p(r,z) = Po{r) exp(—z 2 /2h 2 (r)), valid for a vertically isothermal, hydro- 
static, non self-gravitating disk. We use power-law distributions for the dust surface density 
S(r) = S {r/r )~ p and the scale height h(r) = h (r/r ) 7 where r is the radial coordinate 
in the equatorial plane. The parameter ho is the scale height at a fiducial radius chosen 
to be ro = 100 AU for easier comparison with previous BD and T Tauri disk models. The 
scale height parameters, 7 and h , are free parameters in the modeling and not calculated 
from assumptions of hydrostatic equilibrium. The disk extends from an inner radius r in to 
an outer limit radius r ou t- The central star is represented by a sphere radiating uniformly. 

We consider hom ogeneous spherical grains and we use the dielectric constants of astro- 
silicates (iDraind 120031 ) . The differential grain size distribution is given by dn(a) oc a -3 ' 5 da 



with grain sizes between a min = 0.03wm and a max = 1/mi or 1mm. The mean grain density is 
3.5 g cm -3 . Dust extinction and scattering opacities, scattering phase functions and Mueller 
matrices are calculated using Mie theory. The grain properties are taken to be independent 
of position within the disk. All disk masses discussed below are described as "Total" disk 
mass of gas and dust together. We make no assumptions about the gas distribution in 
the disk, but we assume a ratio of 100 for gas mass relative to dust mass for comparison 
with other disk mass estimates. Of course the SED modeling is only sensitive to the dust 
mass, and even in that case only to dust up to some nominal size of a few times the longest 
observed wavelength. In principle, the two radiative transfer codes should give the same 
result, but because of differences in the exact model grid, there were sometimes modest 
differences. For example, the MCFOST grid in mass was spaced by 0.5 in logio(Mass) from 
10 -6 — 10~ 2 M Q (total disk mass), while the MC3D grid was spaced by 1.0 in logio(Mass) 
from 5 x 10~ 8 - 5 x 10" 4 M . 

As discussed in Paper I, the choice of outer radius for the disks makes essentially no 
difference to the model SED in the spectral range of our study because most of the far-ir 
emission comes from dust at radii inside of 20 AU; so this was left fixed in all the models at 
100 AU. The choice of outer radius also makes only a small difference in the probable values 
of most other disk parameters. This is shown both in Paper I by a comparison of the disk 
parameters f or 2MASS120733 f or outer radii of of 75 versus 25 AU, and by the modeling 
described by iBouv et al.l J2008h for 2MASS04442713. For example, the total disk 



mass is 



driven largely by the total far-infrared/submm flux (see also discussion below and Figure 

m 



Because of the remaining degeneracy in the model-fitting despite our extended SED 
coverage, we present the results largely from the Bayesian probability analysis for various disk 
parameters. In other words, at this point in the model-fitting it is not appropriate to choose 
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a "best-fit" model or model parameter(s) but simply to state the most probable range for 
disk parameters. For example, sometimes the model with the absolutely lowest chi-squared 
implies a parameter value somewhat offset from the peak of its probability distribution, 
since the probability distribution for each parameter is marginalized over all possible values 
of the remaining parameters. This can lead to noticeable differences between some of the 
disk parameters for two models that fit the SED nearly equally well. 

The overall results of the modeling are displayed in several figures. Figures l9"HT2l show 
example model SED's from the set that produces reasonable fits for each object. Figure [T41 
shows histograms of the total of the probability distributions for all the program objects 
as well as the subset detected at 160/zm to show the likely average distribution of these 
parameters. Figure (TS] shows the distribution of the most probable disk masses as a function 
of spectral type for the 43 objects with published spectral types, and finally Figure [TBI shows 
the distribution of probable disk masses versus the 70 and 160/im luminosities for each object. 
We discuss these results below. More refined modeling of individual sources will be pursued 
in a subsequent study where we will report individual disk parameters. 

In Paper I we discussed the comparison between our model results for the first three 
objects and several existing models in the literature for those objects. To our knowledge 
only two additional object s from our s ample, CFHT Tau 9 and OTS 44, have a published 



model. In a recent study iRiaz et al.l (j2012al ) investigated the radial dependence of grain 



composition in brown dwarfs relative to T Tauri stars. Since they used a geometrical model 
of the disk similar to those in our study, it is possible to compare the resulting model 
fits directly. In particular, for CFHT Tau 9 without the additional spectral coverage of 
Herschel, they found an inner disk radius equal to the sublimation radius, a disk inclination 
of order 30-40°, and a flaring index of 1.1 provide a good fit to the overall SE P. These values 



are es sentially identical to the probability peaks in our fits for this object. iLuhman et al. 



(120051 ) have modeled the disk around OTS44, but their model was bas ed on differ e nt disk 



parameters making it difficult to perform a detailed comparison. Finally, IRiaz et al.l (l2012bl ) 
have recently reported a large submm flux in the Herschel SPIRE bands from 2M1207 that 
would imply a m uch larger disk ma ss than found from our modeling, but have since retracted 



that conclusion (IRiaz et al.ll2012d ) 
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6. Discussion 

6.1. Overview of Modeling Results 

Despite the inherent degeneracies in SED modeling without angular size information, 
there are several conclusions that appear relatively firm as illustrated in the figures described 
above. The majority of the disks are best modeled with scale heights h at the fiducial radius 
of 100 AU of order 5-20 AU and radial scale height dependences 7 in the range 1.1 - 1.2. 
Their inner radii are likely close to the dust sublimation radii which are typically of order 
5 x 10 -3 AU for the faintest objects up to ~ 1.5 x 10~ 2 AU for the most luminous. Their 
surface density gradients are best fit with a relatively shallow slope of p ~ 0.5 — 1.0 (X oc r~ p ) 
although this value is not well constrained. There is clearly a very wide range in likely masses 
with a median value of order 3 x 10~ 5 M (0.03 M Jup ), a value which is roughly the same for 
both the objects with spectral types later than M6 and earlier than M6. Unless these masses 
are severely underestimated due to grain growth to sizes well beyond a value to which this 
study is sensitive, this suggests that future giant planet formation would be very difficult in 
these disks. If we assume that spectral type is a rough proxy for brown dwarf mass, then 
Figure IT5l suggests that the ratio of disk mass to sub-stellar mass covers quite a wide range, 
of order a few x 10~ 5 up to ~ 10~ 2 . 

It is well understood that the mass of circumstellar disks is best measured at submm 
or mm wavelengths where the disks are optically thin except in the very central regions 
and the dust is virtually all emitting on the Rayleigh- Jeans part of its SED. The longest 
wavelength observed in our study is 160/im. Although the optical depth in the plane of 
all our model disks is well above unity at 160/aul, Figure [17] shows that for most of the 
objects which are naturally not observed edge-on, the disks are vertically optically thin in 
emission over radial distances where more than half the total flux is emitted. This figure 
also shows that the situation at 70/zm is more problematic with only ~ 1/3 of the total flux 
emitted outside the r = 1 radius. This implies that the 160/im flux may be an approximate 
measure of disk mass, limited mainly by its relative insensitivity to very cool and very large 
dust grains and by the fact that 160/im is not completely into the Rayleigh- Jeans part of 
the dust emission spectrum. Figure [TO] illustrates this possibility by comparing the most 
probable model disk masses to the observed 70 and 160/im luminosities. With the exception 
of the two outliers, both the TWA objects, there does appear to be a reasonably good 
correlation between 160/im flux and likely disk mass, as well as even with 70 /im luminosity. 
No such correlation exists between the total luminosity, for example, and 160/im flux. Even 
though our 160/im data already provide good disk mass estimates, further observations of 
all these objects with ALMA will be extremely important to determine both the disk fluxes 
at (sub)mm wavelengths and disk sizes. 



6.2. Host Cloud Correlations? 



The 46 low-mass objects with disks in our study are located in various star-forming 
clouds as listed in Table [TJ These host clouds have a range of estimated ages as well as other 
less quantifiable differences. For example, the TW Hydra Association is a very low stellar- 
density group, while some of the younger regions have a much higher density of young stars 
and brown dwarfs in their core. In any case, a number of figures in our study that distinguish 
between these host clouds display no clear correlation of any observed or modeled parameter 
with the location of the objects. We have also investigated a simple grouping of objects into 
two bins, those in clouds with estimated ages of order 10 Myr, TWA and Upper Sco, and 
the remainder that are significantly younger. Even with this division which encompasses 
13 objects in the "older" group, we find no significant differences in disk properties. The 
median probable disk mass in the "older" group is logi^Mdisk = —5.0 ± 0.85 M while the 
median for the younger group is logwMdisk = —4.5 ± 1.05 M . 



6.3. Comparison With Young Stellar Objects 



In Paper I we made an initial comparison of the disk properties of our first three program 
objects with young stars that will eventually burn hydrogen, the T Tauri stars. Because 
such stars are brighter than brown dwarfs, there exists a large quantity of observations 
and modeling in the literature to charact erize their disks vi a dust emission, scattering, and 
even gas emission. The recent review by IWilliams &: Ciezal ( 1201 ll ) summarizes the state of 
understanding of disks around T Tauri stars. The disk geometries for such stars are quite 
similar to the geometries found to be likely fits to the disks around our sample of very low 
mass objects within the mutual uncertainties. The typical surface density gradient we found 
of £ oc r~°' 5 to r~ L0 and typical scale height and flaring mentioned above are certainly well 
within the range for models of solar-mass young objects. It is interesting, though, that if the 
disk vertical structure is determined solely by hydrostatic equilibrium with well-mixed gas 
and dust, lower central obj ect masses would le ad to more flared disks around brown dwarfs 



than around T Tauri stars (Walker et al. 



20041 ) which is not obvious in our results. This has 



already been noted by ISziics et al.l (120101 ) in a comparison of median SED's between BD's 
and T Tauri stars. The most significant clear difference between our sample of sub-stellar 
objects and T Tauri stars is the total disk mass inferred from modeling and submm/mm 
observations. Disk masses around T Tauri stars are commonly found to be in the range of 
10~ 3 — 10 _1 M Q , well above our median disk mass of 3 x 10~ 5 M . The most massive disks 
i n our sample are comparable to the masses of a few previously modeled brown dwarfs (e.g. 
( iBouy et al.ll2008l )) and to the low end of the distribution of disks around T Tauri stars. The 
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masses of our sample objects are uncertain, but assuming a median mass ~ few x 1CT 2 M 
for the brown dwarfs, the median ratio of disk mass to central object mass may be a factor 
of 3-10 lower than for T Tauri stars. When the results from the shallow large scale Herschel 
surveys are available, it will be interesting to re-examine this conclusion with the inclusion 
of brighter brown dwarfs. 

7. Summary 

We have used Herschel /P ACS to observe nearly 50 young sub-stellar and nearly sub- 
stellar objects with likely circum "stellar" disks implied by previous Spitzer data, and we 
have increased the number of such objects that are detected beyond 24/xm by a factor of 
~5. We find a wide range in ratios of far-ir to mid-ir flux. Preliminary modeling of the 
disks to fit the observed SED's shows that the dust disk geometries are similar to those of 
their higher mass counterparts, but the range of masses extends well below the masses of 
disks around T Tauri stars, i.e. down below 10~ 6 M . This implies that for most of the 
objects, unless giant planets have already formed, there is not nearly enough mass to form 
them. Our radiative transfer modeling shows that HerscheVs 160/iin band is already long 
enough to provide reasonable disk mass estimates. Although our sample selection precludes 
drawing conclusions about disk frequencies, we have found no clear dependence of disk mass 
or geometry on the age of the star-forming region where the object is located. 
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Table 1: Location of Observed Young Brown Dwarfs 



Cloud 


Number of BD s 


Assumed Distance 


A 1 A 

Assumed Age 


References 


Opmucnus 


3 


125 pc 


1 O TV J 

1-3 Myr 


1 


Taurus 


4 


140 pc 


1-3 Myr 


2 


Chamaeleon I 


12 


160 pc 


1-3 Myr 


3 


Chamaeleon II 


3 


178 pc 


1-3 Myr 


1 


Lup I 


1 


150 pc 


1-3 Myr 


1 


Lup III 


5 


200 pc 


~5 Myr 


4 


Sigma Ori 


7 


360 pc 


~3 Myr 


5 


Upper Sco 


11 


145 pc 


~11 Myr 


6 


TWA 


2 


54 pc 


8-10 Myr 


7 



References. — (1) As adopted by Spitzer c2d, Evans et al. 2007; (2) Kenyon, Gomez & Whitney 2008; (3) 
Luhman 2008; (4) Comeron, F. 2008, 2009; (5) Caballero et al. 2007; (6) Pccaut, Mamajek & Bubar 2012; 
(7) Webb et al. 1999 
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Table 2: Observational Parameters 



Parameter 



Value 



Comments 



OBSID Type 
Wavelengths 
Number of Scan Legs 
Scan Length 
Cross Scan Step 
Scan Angles 
Repetitions 

Peak Intg Time Per Pixel 



PACS Mini-Scan-Map Two Crossed OBSID's 
70/zm, 160/im 
8 
3' 
4" 

70°, 110° Relative to Detector 

7 Per OBSID 

504 sec Per OBSID 



Table 3. Source Summary (Program ID = GTl_pharve01_2) 



Object RA/Dec Center (J2000) OBSID's Obs. Date F„ 70^m F„ 160/«n Sp. Type Refs a 





h m 


s 


< 




i a 










mJy 


mJy 






CFHT-Tau9 


04 


24 


26. 


.5 +26 


49 


51 


1342227059/60 


2011- 


-08- 


-22 


9.5 ± 1.2 


< 6 


M6.2 


14 


KPNO-Tau6 


04 


30 


07. 


.2 +26 


08 


21 


1342227011/2 


2011- 


-08- 


■21 


2.3 ± 1.0 


< 5 


M9.0 


14,15 


KPNO-Tau7 


04 


30 


57. 


.2 +25 


56 


40 


1342227999/8000 


2011- 


-09- 


■04 


3.5 ± 1.0 


< 7 


M8.2 


14 


CFHT-Taul2 


04 


33 


09 


.5 +22 


46 


49 


1342227013/4 


2011- 


-08- 


■21 


1.5 ± 0.7 


< 8 


M6.5 


14,15 


SOri053825 


05 


38 


25. 


A 

.4 


no 
— uz 


42 


41 


1342226721/2 


2011- 


-08- 


•17 


2.0 ± 1.1 


< 7 


M7.0 


7,5,6 


SOri053834 


05 


38 


33 




y 


no 
— uz 


45 


08 


1342226723/4 


2011- 


-08- 


-18 


7.4 ± 0.9 


10.1 ± 5.0 


M4.0 


33,5,6 


SOri053848 


05 


38 


48. 


.z 


no 
— uz 


44 


01 


1342226725/6 


2011- 


-08- 


-18 


2.7 ± 1.0 


< 4 




5,6 


SOri053855 


05 


38 


55. 


A 

.4 


no 
— uz 


41 


21 


1342226727/8 


2011- 


-08- 


■18 


1.9 ± 1.0 
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M5.0 
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SOri053902 


05 


39 


02. 


n 
u 


no 
— uz 


35 


01 


1342228431/2 


2011- 
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-07 


27.9 ± 4.5 


20.0 ± 4.0 


M2-4 


8,5,6 


SOri36 


05 


39 


26 


o 
y 


no 
— uz 


36 


56 


1342228429/30 


2011- 
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-07 


1.5 ± 0.8 


< 5 




5,6 


SOri054004 


05 


40 


04. 


.5 


-02 


36 


42 


1342228363/4 


2011- 


-09- 


-09 


<1.2 


< 3 




5,6 


SSSPM110209 


H 


02 


09. 


.8 


-34 


30 


36 


1342221849/50 


2011- 


-05- 


-29 


9.4 ± 0.6 


8.6 ± 1.8 


M8.5 


30,13,32 


2MJH0703 


11 


07 


03. 


.7 


-77 


24 


31 


1342223480/1 


2011- 


-06- 


-23 


24.6 ± 1.2 


50.1 ± 12.0 


M7.5 


18,13,18,19 


ChaHa9 


11 


07 


19. 


,2 


-77 


32 


52 


1342223482/3 


2011- 


-06- 


-23 


12.2 ± 0.8 
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M5.5 
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11 


08 


04. 


.0 


-77 


30 


38 


1342223484/5 


2011- 


-06- 


-23 


94.4 ± 1.6 


84.1 ± 8.0 




13,18 


IS0138 


11 


08 


19. 


.0 


-77 


30 


41 


1342218699/700 


2011- 


-04- 


-16 


3.4 ± 0.8 


<15 


M6.5 


23,13,18,21 


IS0143 


11 


08 


22. 


.3 


-77 


30 


28 


1342218699/700 


2011- 


-04- 


-16 


13.3 ± 0.9 


10.2 ± 2.0 


M5.0 


23,13,21,22 


ChaHa6 


11 


08 


40. 


.2 


-77 


34 


17 


1342223486/7 


2011- 


-06- 


■34 


14.4 ± 1.0 


<10 


M5.8 


22,13,23 


T37 


11 


08 


50. 


.9 


-76 


25 


14 


1342223468/9 


2011- 


-06- 


•22 


2.9 ± 0.9 


< 3 


M5.2 


23,13,22 


IS0165 


11 


08 


55. 


.0 


-76 


32 


41 


1342223470/1 


2011- 


-06- 


■22 


42.5 ± 0.8 


30.0 ± 10.0 


M5.5 


23,13,18,21 


Cha726 


11 


09 


52. 


.2 


-76 


39 


13 


1342223474/5 


2011- 


-06- 


•23 


10.4 ± 1.3 


<30 


M6.2 


22,13,21,23 


OTS44 


11 


10 


11. 


.4 


-76 


32 


13 


1342223472/3 


2011- 


-06- 


-22 


7.0 ± 1.0 


<10 


M9.5 


24 


IS0252 


11 


10 


41. 


.4 


-77 


20 


48 


1342223478/9 


2011- 


-06- 
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8.7 ± 0.7 


<10 


M6.0 


23,13,18,21 


2MJ11H45 


11 


11 


45. 


.3 


-76 


36 


50 


1342223476/7 


2011- 


-06- 


-23 


3.0 ± 1.0 


<10 


M8.0 


20,13,18,21 


2MJ120733 


12 


07 


33 


.4 


-39 


32 


54 


1342202557/8 


2010- 


-08- 


■10 


7.1 ± 0.7 


< 7 


M8.0 


31,13,27,29 



Table 3 — Continued 



Object RA/Dec Center (J2000) OBSID's Obs. Date F„ 70/mi F„ WOfim Sp. Type Rcfs a 
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.5 
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56 
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.0 
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usdl60958 
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Table 3 — Continued 



Object RA/Dec Center (J2000) OBSID's Obs. Date F„ 70/un F„ 160/xm Sp. Type Refs a 
h m s o i it m jy m jy 



a First reference listed is for spectral type if measured; additional references are for various A < 40/un 
photometry 

References. — (1) Allers et al. 2006; (2) Gully-Santiago, Allers k Jaffe 2012; (3) Bouy et al. 2007; 
(4) Wright et al. 2010 (WISE); (5) Caballero et al. 2007; (6) Hernandez et al. 2007; (7) Rigliaco et 
al. 2011; (8) Caballero et al. 2008; (9) Allen et al. 2007; (10) Comeron 2011; (11) Comeron, Spezzi k 
Lopez Marti 2009; (12) Lopez Marti, Eisloffel k Mundt 2005; (13) Skrutskie et al. (2MASS) 2006; (14) 
Guieu et al. 2007; (15) Rebull et al. 2010; (16) Merin et al. 2008; (17) Alves de Oliveira et al. 2010; 
(18) Luhman et al. 2008; (19) Luhman k Muench 2008; (20) Luhman 2007; (21) Lopez Marti et al. 
2004; (22) Damjanov et al. 2007; (23) Luhman et al. 2004; (24) Luhman et al. 2005; (25) Scholz et al. 
2007; (26) Martin, Delfosse k Guieu 2004; (27) DENIS; (28) Ardila, Martin k Basri 2000; (29) Riaz, 
Gizis k Hmiel 2006; (30) Scholz et al. 2005; (31) Gizis 2002; (32) Luhman et al. 2010; (33) Cody k 
Hillenbrand 2011 
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Fig. 1. — Observed luminosities versus spectral types for the 43 objects with determined 
spectral types in the literature. The symbols are color coded so that the two star-forming 
regions that are significantly older than the others, U Sco and TWA, are in red. 
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Fig. 2. — Histogram of observed luminosities (1-200/xm) for all 47 objects and histogram of 
spectral types for the 43 objects with determined spectral types. 




Fig. 3. — Images of first 10 objects in Table |3j Dotted circle is 30" in diameter. 




Fig. 5. — Images of next 10 objects in Table [3) 




Fig. 6. — Images of next 10 objects in Table [3j 
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Fig. 7. — Images of last 6 objects in Table [3j 
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Fig. 8. — Overplotted SED's normalized to 1.6/im. Objects with S/N > 2 at both 70 and 
160;um are shown with fluxes plotted as red circles and red SED's. Objects with S/N > 2 
only at 70/mi are shown plotted in green, and remaining objects without reliable detections 
at either wavelength are in blue. Plotted fluxes without circles are only upper limits or low 
S/N. The solid black curve shows a typical SED of a bare BD photosphere (SSSPM110209 
from Paper I). The possible Class I object, Chall0804, is the solid blue line with the highest 
values of normalized 70 and 160//m fluxes; note flux values for his object are upper limits 
shortward of 3.6/im as measured from the 2MASS survey images. 




Fig. 9. — SED's of the first 12 Class II objects in Table [3j Representative model fits are 
shown with dashed lines, and Spitzer IRS data are shown in red. A typical bare photosphere 
is shown in solid, normalized to the observed 1.6/im flux. 
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. 10. — SED's of the second 12 Class II objects in Table [31 
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Fig. 11.— SED's of the third 12 Class II objects in Table H 
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Fig. 12.— SED's of the last 10 Class II objects in Table EJ 
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Fig. 13. — Color-color diagram showing far-ir color {vF u 24/70) versus mid-ir color [yF v 
3.6/24) as a function of star-forming cloud in which the object is found. There is no 
clear trend with region, suggesting that there is no obvious dependence of outer disk struc- 
ture/mass versus age. 
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Fig. 14. — Histograms of the total of all probability distributions for the disk parameters 
illustrating the likely range of these parameters in the total sample. The solid line is for the 
entire sample and the dashed line is for just the objects detected at 160/im. 
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Fig. 15. — Most probable disk mass versus spectral type for the 43 objects with published 
spectral types. 
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Fig. 16. — Comparison of the 70 and 160/xm luminosities with the model-derived most 
probable disk masses illustrating the utility of the longest PACS wavelength for estimating 
disk masses. 




Fig. 17. — Upper panel - plot of cumulative emitted flux as a function of radius at several 
wavelengths for a typical disk model with M disk = 2 x 10~ 4 M , R in = 0.05 AU, a max = 1mm, 
i = 30°, and geometry equal to the likely values in Figure HH Lower panel - optical depth 
perpendicular to the disk as a function of radius at the same wavelengths, illustrating that 
at 160/im and longer, the disks are typically optically thin over most of the radii emitting 
at those wavelengths. 
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